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We report on the progress of the pentaquark searches by the STAR collaboration in p+p, d+Au, 
and Au+Au at snn — v / 200 GeV collisions through one of the decay modes of + — > p + K®. 6> + 
state is an exotic baryon with strangeness S = 1 and is the lightest isospin member of the expected 
antidecuplet. We compare our techniques for the pentaquark search with those for the short-lived 
resonances. These results were presented as a poster at Quark Matter 2004. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

The first observations of the 9 + pentaquark, a five 
quark bound system of uudds, have been reported in 
photon- nucleus and kaon- nucleus reactions Q, 0, . The 
presence of this state was predicted by R. Jaffe with mul- 
tiquark bag models and later by D. Diakonov et al. 
using chiral soliton models of baryons [(| . 

The high energies and particle densities resulting from 
collisions at the Relativistic Heavy Ion Collider (RHIC) 
are expected to favor pentaquark production. The data 
that we analyze were taken by the STAR (Solenoidal 
Tracker At RHIC) experiment, one of the four relativis- 
tic heavy ion experiments at RHIC. The large acceptance 
of STAR'S Time Projection Chamber (TPC) is ideal for 
such rare particle searches. The short lifetimes pre- 
dicted for pentaquarks require that a mixing technique 
be used to reconstruct pentaquarks via their decay prod- 
ucts. This technique has already been used successfully 
by STAR to reconstruct and study resonances 0, @ • 



II. ANALYSIS AND PARTICLE 
IDENTIFICATION 

Charged daughter particles are identified by the mo- 
mentum and the energy lost per unit length, dE/dx, 
measured with the Time Projection Chamber of STAR. 
For long lived particles (ct ~ few cm) such as K®, A 
and 5, the decay vertex topology information is used 
for their identification. This method cannot be used for 
pentaquarks since they decay strongly (ctq ~ fm). An 
alternative method, a mixing technique, is used success- 
fully to identify short lived particles such as resonances 
( CT E(i385) = 5 fm) and can be used to search for pen- 
taquarks and dibaryons. 



In the mixing technique, for example £(1385) — * A + 7r, 
we identify first the A's by their decay vertex topology, 
due to their long lifetime (ct\ — 7.89 cm). This A candi- 
date is then combined with a it to get £(1385). The back- 
ground is described by combining 7r's from one event with 
the A's from another event. Similarly for pentaquarks 
such as + — > K® + p, we first identify the K® via its 
decay vertex topology and then combine this K® candi- 
date with a proton candidate in the same event to extract 
the signal and in a different event to describe the mixing 
background. 



III. STATUS OF CURRENT STUDIES 

A. Monte Carlo Studies 

To study the decay mechanism and optimize the ap- 
plied cuts, we use Monte Carlo simulations. In this study, 
one Monte Carlo + pentaquark is chosen from a ther- 
mal exponential distribution with T invs i ope = 250 MeV 
for a rapidity | y |< 1.5 and is embedded in a single real 
p+p event after a full TPC simulation. The chosen in- 
put width, 10 MeV/c 2 , and the input mass, 1.54 GeV/c 2 , 
are consistent with the observed mass and width of + 
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FIG. 1: Invariant mass spectrum of the simulated Monte 
Carlo input for the + and pr vs rapidity of the <9 + that 
can be identified with the TPC. 
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HHI3. In Fig.HJthe invariant mass spectrum of simu- 
lated 8 + and the TPC acceptance for the corresponding 
simulated particles can be found. In Fig. the invariant 
mass spectrum of the reconstructed + are presented. 
We find that ~ 3% of these Monte Carlo + 's are suc- 
cessfully reconstructed with this technique. The recon- 
structed width and the mass is consistent with the Monte 
Carlo input. We can study the decay properties with the 
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FIG. 2: Invariant mass spectrum of the Monte Carlo simu- 
lated Q + embedded in real p+p events with the mixing tech- 
nique. Black solid line is the signal and red dashed line is the 
mixed event background. The simulated signal can be clearly 
seen. 

simulated tracks such as the momentum distribution of 
the decay daughters. Fig. [3|and Fig. 0]are the momen- 
tum distributions of the K® and p respectively. With 
this study we can apply momentum cuts to increase the 
signal-to-background ratio. This ratio can be increased 
clearly by selecting protons with momentum less than 
1 GeV/c. Further detailed studies on other variables are 
needed to optimize the signal-to-background ratio. 
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FIG. 4: Proton momentum distribution: Black solid line is of 
the accepted p after the dE/dx cut and red dashed line is of 
the decay daughters of the Monte Carlo simulated + for the 
same number of events. 



In these figures the invariant mass spectra of £(1385) + , 
£(1385)~ and their antiparticles are summed to improve 
the statistics. Similarly Fig. is for the Au+Au central 
collisions. For this figure only £(1385)+ and £(1385)" 
are added. 
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FIG. 5: Background subtracted invariant mass spectrum of 
£(1385) in p+p collisions. 
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FIG. 3: Kg momentum distribution: Black solid line is of the 
accepted K° after the topological cuts and red dashed line is 
of the decay daughters of the Monte Carlo simulated <9 + for 
the same number of events. 



3 10000 



5000 




1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 

m Inv ( A+n ) [GeV/c 2 ] 

FIG. 6: Background subtracted invariant mass spectrum of 
£(1385) in d+Au collisions. 



B. Mixing Technique Works 

£(1385) —>■ A + ^bachelor is identified in p+p, d+Au, 
and Au+Au collisions with the mixing techniques. Fig[5] 
and Fig are the preliminary, background subtracted 
invariant mass spectra of the 3 — > A + tt and £(1385) — > 
A + 7r for p+p and d+Au collisions at ^Snn = 200 GeV. 



The 3 peaks are fit with gaussian distributions and the 
£(1385) peaks are fit with a Breit-Wigner distribution. 
Antiparticle-to-particle ratios of 3 and £(1385) are ob- 
served as 0.89 + 0.04 and 0.90+0.07 respectively in p+p 
collisions. (In addition to the given statistical errors, 
these ratios also have a 20% systematic error due to the 
mixed event background normalization.) These values 
are consistent with what we observe for other hadronic 
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FIG. 7: Background subtracted invariant mass spectrum of 
£(1385) in Au+Au central collisions. 

particle ratios at RHIC in p+p collisions, such as the an- 
tiproton to proton ratio which is p/p ~ 0.8. We expect a 
similar production of pentaquarks and their antiparticles 
at RHIC. A very preliminary invariant mass spectrum 
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FIG. 8: Invariant mass spectrum of the Q + in p+p collisions. 
Black solid line is the signal and red dashed line is the mixing 
background for the upper figure. The lower figure is the mixed 
event background subtracted invariant mass spectrum. The 
expected position of the mass of G is at 1.54 GeV/c 2 . 

of the © + in p+p collisions is shown Figure 00 Further 
studies are needed to optimize the signal-to-background 
as no strong signal is observed. The invariant mass spec- 
trum that we observe for the + in p+p collision events is 
consistent with our predictions for the significance of the 
signal given our current statistics. Details can be found 
in flIII CI The signal-to-background depends highly on 



the selection of events and applied cuts. To improve cuts 
and understand the decay mechanism we can do more 
detailed simulation studies. For Au+Au and d+Au col- 
lisions similar invariant mass spectra are observed. 



C. Feasibility Studies 

Assuming the + production is 10 — 100% of the 
A(1520) in p+p collisions, one can predict the yield of 
+ . Preliminary dN/dy of A(1520) at mid-rapidity is 
0.004 per event in p+p collisions 0,0]. There are 8 
Million p+p events available for this analysis. This cor- 
responds to a production of ~ 30000 A(1520) and a pro- 
duction range of ~ 3000 to ~ 30000 + 's in these p+p 
events. As the efficiency of the mixing technique is ~ 3% 
and the branching ratio of the — > K% + p is ~ 25% 
(assuming that the branching ratios of — > K° + p and 
—>■ K+ + N are 50% each), 20-200 of the 0's can be 
found with the mixing technique. The background pairs 
per event in the 1.54+5 MeV mass range is 3200. This 
corresponds to a significance of 0.25 to 3 for the signifi- 
cance defined as ,„ 5 Sl 9 nal _ Similarly one can 

V 2x Background+bignal ^ 

repeat the same study for Au+Au and d+Au collisions 
and correspondingly predict a significance of 2-7 for 1.5 
Million Au+Au events and 1-16 for 10 Million d+Au for 
the predicted production of one + per rapidity per col- 
lision 0, 0, [n|. To estimate the yield for the d+Au 
collisions we assume N part scaling. The mean number of 
participants in d+Au is 8, in p+p it is 2, and in Au+Au 
it is 350 for the most central collisions. The lower limit is 
obtained from p+p scaling while the upper limit is from 
Au+Au yield estimates. 



IV. CONCLUSION 

Preliminary acceptance and efficiency studies show 
that we should be able to find pentaquarks at the few % 
level. Resonances can be clearly reconstructed via event 
mixing techniques in p+p, d+Au and Au+Au central 
collisions. Optimization of cuts to improve the signal 
over background is in progress. There is a possibility of 
measuring the anti-pentaquarks at RHIC since the an- 
tibaryon to baryon ratio indicates a nearly net baryon 
free region [1J|. 100 Million events for Au+Au collisions 
at y/sNN = 200 GeV are planned to be collected in 2004, 
Run 4, which has just started. This would be 70 times 
the currently available data. The predicted significance 
of identifying + is at a range of 20 to 84 with the forth- 
coming data. We expect to put a more definite upper 
limit to the yields and production mechanisms of the pen- 
taquarks in Au+Au collisions with the 2004 run which 
will finish collecting data at the end of May 2004. 
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